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The addition of gauge singlet fermions to tiie Standard Model Lagrangian renders the neutrinos 
massive and allows one to explain all that is experimentally known about neutrino masses and lepton 
mixing. At the same time, the gauge singlet fermion decays in the early universe produce a lepton 
asymmetry, which is converted to a baryon asymmetry via Spharelon processes (leptogenesis). On 
the other hand, the addition of a gauge singlet scalar to the Standard Model yields a thermal dark 
matter candidate through interactions between the Higgs boson and the gauge singlet scalar. By 
imposing a Z2 symmetry on the gauge singlet scalar and one of the gauge singlet fermions, we can 
have viable dark matter candidates and new interactions coupling the Z2-odd scalar to the Z2-odd 
fermion, which can lower the leptogenesis scale (and the reheating temperature) to 0{TeV). 

I. INTRODUCTION 

Gauge singlet fields are a simple but very interesting form of physics beyond the Standard Model (SM). A gauge 
singlet scalar(S') can be a thermal dark matter candidate by having S couple to the Higgs boson (for example, see 
In order to have it be stable or long-lived, a symmetry might be imposed or any decay channel has to be suppressed 
by a high energy scale, usually comparable to the Planck scale. 

On the other hand, gaugc-singlct fermions, known as right-handed (RH) neutrinos, can explain the observed tiny 
neutrino masses Q via the type-I seesaw mechanism Q. RH neutrinos can also accommodate the observed baryon 
asymmetry [3| via thermal leptogenesis Q. The mechanism of leptogenesis satisfies the three Sakharov's conditions 
Q (i) baryon number violation, (ii) C and CP violation, (iii) deviation from thermal equilibrium. A lepton number 
asymmetry is generated via the decay of heavy RH Majorana neutrinos, which is converted into the baryon number 
asymmetry through Spharelon processes . In order to be the main source of the baryon asymmetry, the mass scale of 
the heavy RH neutrinos must typically be larger than 10^ GeV [8|[, which requires a high reheating temperature. Such 
a high reheating temperature leads to gravitino overproduction. [9| in the context of supersymmetry. There are many 
ways to avoid gravitino overproduction. Resonant leptogenesis [lO|, for instance, assumes the limit rriN^ — mNi <C m^Va 
so that rriNi and m^Va can be as low as of order TcV. There is enhancement to the lepton asymmetry by taking into 
account flavor effects [ll| . In [l2j , O , and references therein, there are more detailed discussions on the solutions to 
the gravitino overproduction problem. 

It is intriguing to combine these different ideas together, i.e., to have the gauge singlet scalar and fermions at once 
in the theory. This is done, e.g., in [T^], which aims at lifting the tension between X-ray bounds and the Lyman-a 
bounds in the Dodclson-Widrow (DW) model [l^, where the dark matter is the RH neutrino that is generated from 
neutrino oscillations. The constraint from X-rays puts an upper bound on the mass of the RH neutrino since it can 
decay into an active neutrino and a photon via loop diagrams and the decay rate is proportional to the mass of the RH 
neutrino. At the same time, Lyman-a puts a lower bound on the thermal-averaged momentum of the RH neutrino 
which equivalently is a lower bound on its mass. For more details, see In the context of the DW model, the two 
constraints can not be satisfied at once. To reduce the tension between the two bounds, has both gauge-singlet 
scalar S and gauge-singlet fermions N, where S couples to N. S decaying into N provides an extra production 
mechanism of N in addition to neutrino oscillations. A^'s from the decay of S get red-shifted when the universe 
cools down; therefore, they have a lower thermal average momentum than those generated by oscillations. Hence 
the Lyman-a bounds become weaker. In [l6j . a similar setup is employed. With the help of an unbroken Abelian 
family symmetry, Gtamiiy, under which the full Lagrangian is invariant, all particles carry the charge of Gfamiiy The 
gauge singlet scalar, S, couples to gauge singlet fermions, A^'s, with a nontrivial structure and, at the same time, the 
vacuum expectation value of S provides masses to A^'s, in addition to the Majorana mass terms, in such a way that the 
Maki-Nakagawa-Sakata (MNS) matrix (or the Pontecorvo- Maki-Nakagawa-Sakata (PMNS, or MNSP) matrix) yields 
^13 = and some regions of the parameter space with the weak washout effect, which can yield thermal leptogenesis. 

In this paper, we manage to add more structure into theory, i.e., an extra discrete symmetry, so that the framework 
can provide dark matter candidates, alleviate the gravitino problem of leptogenesis and make a connection between 
the baryon asymmetry and dark matter. To be more specific, we introduce a gauge-singlet scalar S and several 
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gaugc-singlct fermions N and impose a Z2 symmetry on S and one of N's. In this situation, the hghter particle of S 
and N charged under the Z2 symmetry can be the dark matter and the Z2 symmetry guarantees the stability of the 
dark matter candidate. For leptogenesis, loop-diagrams with gauge singlets running inside give rise to the required 
strong and weak CP phases for the generation of the lepton asymmetry. We found that the contributions to the 
lepton asymmetry from the new interactions can be a dominant source of matter-antimatter asymmetry without the 
problem of gravitino overproduction. [l7| uses a similar concept with the type-II seesaw mechanism. 

This paper is organized as follows. In Sec|lTl we describe the formalism and the particle content. In Sec lIIIl a 
detailed analysis of leptogenesis from new interactions has been displayed. The discussion on dark matter is in Sec lIVI 
and we conclude in Sec. |Vl 

II. THE FORMALISM 

The full Lagrangian can be written in the following way, 

C^CsM + Ci+C2+C3+Ci, (III) 

where 







(11.2) 


C2 




(11.3) 


C3 


= -V{H, S) ml\H\' - Xh\H\' - \mlS' - \xsS' - ^\H\'S^ 


(11.4) 


u 


= -XiSN'i^^—^N + h.c. 


(11.5) 



CsM refers to the SM Lagrangian and 5 is a real singlet, whose vacuum expectation value (VEV) is zero, and its mass 
to| = ml + Xhs{H)^- We introduce RH Majorana neutrinos, iV,, which generate tiny masses for active neutrinos, 
and N ^ which provides a new mechanism for leptogenesis. Here, we are interested in the situation with 3 active, 3 
RH neutrinos, Ni, and one extra gauge-singlct Majorana fcrmion, N . The Lagrangian preserves a discrete symmetry 
under which S — >■ —S and TV — > —TV and the other particles remain unchanged. Therefore, there are no cubic terms 
for S and no yHLaN, i.e., is not responsible for the see-saw mechanism. 

III. LEPTOGENESIS 

A. Baryon and lepton number 

We start by discussing how large the lepton asymmetry, e, should be in order to explain the observed baryon 
number asymmetry. In [l^ jTsj . the relation between baryon (B) and lepton number (L) before and after sphaleron 
conversion is 

B{tf) ^ ciB - L) = -cLiU), (111.6) 

where subscripts, i and /, refer to before and after the sphaleron process, respectively. Cs — (SA^/ + ANh)/ {'22Nf + 
13Nh)* in which Nf is the number of generations, Nh is the the number of Higgs doublets and B — L is anomaly-free^^ 
so that we do not need to assign a time index to it. From the Wilkinson Microwave Anisotropy Probe (WMAP) data 

r]B = — = moW-^" = 273.9 Qbh^ 10"^" = (6.23 ± 0.17) x lO-^°. (111.7) 



* In the SM, = 28/79 and = 8/23 for the MSSM. 

t In the context of the SM with additional RH neutrinos, the number of RH neutrinos for anomaly cancelation has to be three if _B — L 
is gauged. In this paper, B — L is a global symmetry and there is no constraint on the number of RH neutrinos. 
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FIG. 1: Contribution to leptogenesis from new interactions. 




A similar quantity to 77 is y, which is the number density in a eomoving volume, defined as 

y = (111.8) 

s 

where n is the number density and s is the entropy density. Obviously, the relation between 77 and Y is, 

Yb=Vb^. (III.9) 
s 

Assuming the entropy and baryon number per eomoving volume remain constant from tf to the time of BBN, ^bbn 

and to that of matter-photon decoupling, tvvMAP, then YB(tf) = YB(tBBN) = ^bI^wmap)- Given that at iwMAP, 
s = 7.04 n-y, we have 

^^(i^) = r^B^ = 8.8 X 10"", (III.IO) 

and 

Yl{U) = --Ysitf) = -2.53 X 10-1°, (III.ll) 
where we assume Nh = 2 and Nf = 3 for the MSSM. 



B. leptogenesis with new interactions 



From the new interactions shown in Fig. [U and the original seesaw ones in Fig. [51 we have 



(III.12) 



4 



where mjV; is the mass of Ni (i = 1, 2, 3) and, from |20| . 

f{x) = (\—^ - (1 + x)hi(ii^)) ^ _^(when X » 1) (III.13) 
\ 1 — X ^ J ^ 

and we do not consider the degenerate situation, in which toat^- — m^r. ^ rriNj 

To simphfy the expressions, we use the Casas- Ibarra parametrization [2lj, which separates the high energy physics 
from the low energy one, and fohow the procedure from Q, where they have obtained the bound of rriN-^ > 10^ GeV 
for the standard typc-I seesaw. First, some quantities are defined as follows, 

^mNi 
^?7ijv3 



D 



RR' 

Via 




ijJKD^RD^U^^, (III.14) 



where rrii,. are the active neutrino masses with mjy, > m^^ > m^^ and U is the MNS matrix. 

For pedagogical purposes, we repeat the derivation of the Davidson-Ibarra bound on ei with corrected coefficients 
for the original leptogcnesis Q . From the second term in Eq. ()III.12| , wc obtain 




^yim(yyt)?,^ 
1287r2 ^ '^' mti- 

1^1 



128^2^""^'"' 



^ jY,ml^mlluv{Rl^), (III.15) 



1287r2(if 

where we have used Eq. (1111.14^ and the fact that {yy^ )ii is real. We have, 

_ r(jVi ^ H+1-) - r{N, ^ H~l+) ^ 3 m^. Eg rnljmjRlJ 

T{Ni^H+l-)+T{N,^H-l+) 167: (H)^ Ea"^-J^iaP ' ^ ' ' 

For simplicity, we assume that light neutrino masses are hierarchical, i.e., m^.^ ^ ™i^2J '^i/u 

3 rriN, mlJlmjRl^)] 
ei < TH~7m2 \u \2 ^ (,111. 17j 

where equality holds when |Im(_Ri3)| ^ |Re(i?i3)|. Thus, wc obtain 

mjvi > i^^^^leil ~ lO^GeV, (III.19) 
3toi/, 



■f Even without mass degeneracy, the self-energy diagram, (b) in Fig. [2] contributes to the Icpton asymmetry. See [20II for more detail. 
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where we use the fact that the observed baryon asymmetry is due to leptogenesis to infer |ei|. We will discuss this in 
more detail later. 

Now, we turn to the new contribution to ei from new interactions, the first term in Eq. (jIII.12p . Ai are assumed 
real for simplicity, and 



J2 MyLy^c.) - E - — Im^ " ^HRIM- (ni.20) 



Therefore, the total Icpton number asymmetry from the decay of Ni is, 
r(A^i ^ H+1-) - T{Ni H-1+) 



V{Ni H+1-) + T{Ni H-1+) + T{Ni NS) 



(III.21) 



1 / Ai A, m^y ^mjv, tun, 3 ^ 2 2 . f r,2 ^ 



where we use Eq. ()III.12P , (jllLlSp . andJIIL20|) . 

Note that Yl = tr{Y'^^ (T ^ toati) [22], where 77' is an efficiency factor which measures the wash-out effect and 
Y^'(T ^ niNi) — 135C(3)/(47r'^ q»), where is the number of relativistic degrees of freedom in thermal equilibrium, 
g^, ~ 230 for the MSSM. From [l3|, 77' ~ m^,/maa in the strong wash out scenario when m > m^, and 771^0, > m*, 

where - 10"^ eV and 77i ee ^„ Tfi^a = Ea 87r|^r(iVi ^ H+l-,H'l+) = ^(A? + (yy^)ii). To sum, we have 
135C(3) m.mjvi (^^ Vet" V"'".™«i™"= Im(i?,„i?*„) + ^^-3^ ,7i^^TO2^Im(i?2^)^ 

where we approximate iiiaa by rh. 

If Xi = 0, ry' = 1, 5* = 230, Cs = 8/23, and = y/Am^ in the MSSM and requiring the observed baryon 
asymmetry coming from leptogenesis, we have 

= < 135C(3)^m^ 

^ ^ - 4TT^g* 167r (iJ)^ ' ^ ' 

which infers ?71n^ > 10^ GeV, that is consistent with the result in 0]. 

From Eq. piL22p . it is obvious that m^Vi could be smaller than 10^ GeV if one were to increase A^ or Ria-^ From 
now on, we focus on the case with only two RH neutrinos, A'^i and N2, for simplicity. The generalization to more RH 
neutrinos is straightforward. 

Increasing A2(^ (2/2/^)22) will cause N2 to depart from thermal equilibrium at later times and the relic density of N2 
is roughly exp{— 771^2 /To) *n-y at decoupling, where Td is the decoupling temperature of N2 and is fully determined by 
A2. However, as long as mjvj 3> rriNi ^^d A2 is large enough, the relic density of N2 is too low to have any impact on 
the lepton asymmetry created by m^^ at T ~ rrimi- For instance, we choose m^^ ~ lOmN^^, Xf ~ (yy^)ii ^ (h)^^ 
and = 230, 

8.8x10"" 135C(3) m^TTiAT^ 1 A2 mj^^ y/W^T^nN^rn^^ 



135C(3) m^mNi {H^ A2 rn%^ ^mN^mN^m^^ 



Then, it is easy to show. 



S Increasing Ai or Ric, would not allow one to reduce mjvj , since they will be canceled by the denominator in Eq. (IIII.22I I. 
^ This choice implies the strong wash-out scenario. 
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In a dramatic situation, we can have rriN^ around O(lOTeV) by having Ai ^ 10~^ and A2 ^ 1. 

On the other hand, increasing i?2a yields larger e without changing 77', which is independent of R2a- However, 
the estimate rj' ^ m^./maa, takes into account the effect of inverse decay (/ + H — > iVi) only. By increasing i?2a, 
AL = 2 {l^H^ <-> l^H~ and /+/~ O- H~ H^) interactions will become important as well. To be more precise, from 
Eq. (|III.22p . for MSSM in the limit of 3> tumi ^■iid strong wash out region, ~ > {yy'^)ii > (h)"^ ' 
constraint on tojv^ is 

For example, if we would like to make m^r = 10^ GeV, Im(i?2a) has to be 10^^ ~ 10"^ for m^v, = lOrriN. 

and ^ 1. The ratio of AL = 2 interactions mediated by N2 to those mediated by A^i will be roughly 

l^2a/^iaP ^ 10^^, which is extremely large and has to be carefully considered in the calculation of In other 
words, increasing i?2a may not be an efficient and applicable way to lower mjvj. 



IV. THE DARK MATTER 



As mentioned before, due to the discrete Z2 symmetry, the lightest of S and N can be a thermal relic. We first 
calculate the relic abundance for each of them, respectively, and then discuss if any of them can be the dark matter 
and at the same time the low-scale leptogenesis is viable. 



A. S as dark matter(mjv > ms) 



From [23| , we know the relic density of the dark matter is determined by the annihilation rate into SM particles in 
thermal equilibrium at the time of decoupling. A rule of thumb is that it decouples when the interaction rates with 
SM particles are smaller than the expansion rate of the universe. To be more quantitatively precise, the Boltzmann 
equation should be used (sec chapter 5 in [1^). With the help of [l|, where a complex scalar is assumed, wc can 
calculate {av) for SS HH , SS ^ W+W- , SS ZZ and SS f+f', respectively. 

•SS HH 



• SS w+w- 



>SS^ ZZ 



iav 



1287rm| \ m|, 



1/2 



IGnmliiAml-mD^+mlrl) 



m 
m 



1/2 



1/2 



(IV.26) 



(IV.27) 



(IV.28) 



• SS^f+f~ 



3/2 



ng^ ((4m| - m2)2 + ^2^2) 



(IV.29) 



where we have an extra factor 1/2 compared to [l[ due to the fact that S is a real scalar field in our case. The relic 
density of S is given by (see Eq. (2.7) in [l|). 



PS giTy) 



K 



7 ^ 



3xf3 
2 



Pc g{Tfs) TjXfsiav) pc 1 



(IV.30) 



where ps is the energy density of S, pc is the critical energy density of the universe, T/s is the freeze-out temperature 
for S^ Xfs = Tfs/ms, is the present photon temperature, g{T) is the the number of relativistic degrees of freedom 
around temperature T, and K ~ {4t:^ g{T) /AbM^i)^/^ . An extra factor 1/2 is again because of S being a real scalar. 
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FIG. 3: Interactions determine the relic abundance of A'^. 



N as dark matter(mg > mjv) 



The main interactions responsible for the rehc density of N are as shown in Fig. [3] Keep in mind that only N and 
S are odd under the discrete symmetry. 

It is expected that (a) in Fig. [3] is the dominant contribution although it is loop-suppressed since (6) is of order 
which is small because we focus on the situation m^Vi < 100 TcV, and there is a phase suppression due to the 
four-body decay. Therefore, we have 



Vai I 



X4 \2 



T 



20487r5 



''Ni 



TUN 



1- ) ( / dxlog 



'TV 



1/2 



''Ni 



(IV.31) 



where A = rrig + x'^E'^^ — ^^cm we have assumed mpfi 3> ms > niiq . The relic abundance of N is given by Eq. 
(|IV.30|) with S replaced by N. 

In principle, if ms ^ ttln, we have to consider co-annihilation interactions, i.e., S + N ^ Ni ^ + , whose 
amplitude squared is proportional to y^j. It is small compared to other annihilation channels because, again, we are 
interested in the situation of mjV; < 100 TeV. 

In summary, for S being the dark matter, (av) is mostly determined by Xhs, which is a free parameter from the 
point of view of Icptogencsis while for iV, (av) is determined by Xi for Ni running in the loop. On the other hand, 
with the help of the large A2 {N2 propagating inside the loop), we can have the correct abundance for N. Fig. |3] shows 
the allowed regions of A2 and Xhs for generating the right dark matter density and having low-scale Icptogencsis. In 
this situation, we have to push both Xi and Xhs toward the strongly-coupled region. 
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FIG. 4: The blue band represents the region of the parameter space of A2 and Xhs, which gives the right thermal relic abundance 
for A'^i and has successful low-scale leptogenesis, where {rriH, m-iv, ms, ttini , m-jva) ~ (100, 150, 200, 3000, 15000) in units of GeV 

V. CONCLUSIONS 

In this papaer, we propose a simple and economical model, which can accommodate both leptogenesis and the dark 
matter, by introducing the new scalar S, which couples to the Higgs boson and RH Majorana neutrinos (Ni), and N. 
We impose a Z2 symmetry under which S and N are odd and the rest is even; therefore, the lighter of S and N could 
be the dark matter. By increasing the coupling, A2, and having 7V2 much heavier than A'^i, we can easily increase the 
efficiency of generation of the lepton asymmetry without having a high reheating temperature and A^i can be as low 
as O(TeV). Note that all of the estimates of the lepton asymmetry are based on the one-flavor approximation instead 
of three flavors (e, ^ and r). Taking into account flavor effects, the efficiency factor, 77', will be modified. From [l^, rj' 
is enhanced by one or two orders of magnitude or remains the same order of magnitude compared to that of one-flavor 
approximation, which implies the estimate done before remains valid. 

Finally, there have been studies of constraints on S as the thermal relic via the decay of the Higgs Boson into S, if 
kinematics allows, or the elastic scattering between S and nuclei, see for example, [2^ . 

As for N as the dark matter, it can be produced in pairs via the Higgs boson and S. If the produced S is on-shell, 

the decay width of S into and a light neutrino is of order (^TT^) ^ 10~^^ GeV, which roughly corresponds to 

the decay time 10"^'^ sec, which implies it can happen inside a detector. However, it is very challenging to identify 
that a new state (other than S) has been produced since both the N and the light neutrino would escape the detector. 

At the time of writing, we notice that, in [25j . they have demonstrated how one may obtain leptogenesis and 
neutrino mass generation in the see-saw picture without invoking a new mass scale far beyond that of electroweak 
symmetry breaking by introducing more than one Higgs doublet family. 
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